Spectral functions of two-band spinless fermion and single-band spin-1/2 fermion 

models 

E. Orignac 
Laboratoire de Physique de I'ENS-Lyon, CNRS UMR 5672, 46 Allee d'ltalie, 69364 Lyon Cedex 07, France 



O 

(N 

D 

00 
(N 



M. Tsuchiizu and Y. Suzumura 
Department of Physics, Nagoya University, Nagoya 464-8602, 
(Dated: September 29, 2011) 



Japan 



We examine zero-temperature one-particle spectral functions for the one-dimensional two-band 
spinless fermions with different velocities and general forward-scattering interactions. By using the 
bosonization technique and diagonalizing the model to two Tomonaga-Luttinger-liquid Hamiltoni- 
ans, we obtain general expressions for the spectral functions which are given in terms of the Appell 
hypergeometric functions. For the case of identical two-band fermions, corresponding to the SU(2) 
symmetric spin-1/2 fermions with repulsive interactions, the spectral functions can be expressed in 
terms of the Gauss hypergeometric function and are shown to recover the double-peak structure 
suggesting the well-known "spin-charge" separation. By tuning the difference in velocities for the 
two-band fermions, we clarify the crossover in spectral functions from the "spin-charge" separation 
to the decoupled fermions. We discuss the relevance of our results to the spin-1/2 Hubbard model 
under a magnetic field which can be mapped onto two-band spinless fermions. 






I 

o 
o 



en 
> 

m 

o 



% 



PACS numbers: 71.10.Pm, 71.10.Fd, 79.60.-i, 67.85.Lm 



I. INTRODUCTION 

The physical properties of interacting fermions in three 
dimensions can be described using the Landau Fermi liq- 
uid concept of fermionic quasiparticles with renormal- 
ized masses and weak effective interactions ii^ In one 
dimension, the Landau Fermi liquid concept is not ap- 
plicable, and interacting spin-1/2 fermions are described 
within the Tomonaga-Luttinger (TL) liquid concept, in 
which quasiparticles are replaced by collective density (or 
charge for electrons) and spin excitations propagating in- 
dependently with respective velocities Up and Ua, the so- 
called spin-charge separation phenomenoni^^— Moreover, 
in the TL liquid, the single-fermion excitations are unsta- 
ble and decay into collective modes, leading to the sharp 
suppression in the density of states.^ Originally, candi- 
date systems for the observation of TL liquid physics 
have been quasi-one-dimensional conductors such as the 
first molecular conductor TTF-TCNQ,^'!" the Fabre and 
Bechgaard salts ^i the blue and purple bronzes, and 
carbon nanotubes.^^ Angle-resolved photoemission spec- 
troscopy (ARPES) experiments in quasi-one-dimensional 
conductors have been used to measure electron spectral 
functions and probe their TL liquid features i^ii^"— Self- 
assembled one-dimensional metallic chains of atoms on 
semiconductor surfaces have also been considered as po- 
tential candidates for the observation of TL liquid and 
spin-charge separationj^i^ Initially, the ARPES mea- 
surement for gold atom chains on silicon surfaces re- 
vealed a double band structure and the spin-charge sep- 
aration was suggested for its origiui ^^'^^ However, recent 
spin-resolved ARPES experiments on this system^- have 
shown that the double band nature originates in the spin 
splitting caused by the Rashba effect. Besides ARPES, 
magnetotunneling measurements between a wire and a 



two-dimensional electron ga o^^i^^ or between two wires^ 
are also sensitive to the TL liquid features of the spec- 
tral functions. Experiments on quantum wire o^^'^^ have 
partially confirmed the presence of TL liquid effects in 
tunneling current measurements. 

More recently, it has been proposed that ultracold 
atomic gases were also candidates for the observa- 
tion of the TL liquid and spin-charge separatioui^"— 
Indeed, atom trapping technology, either optical^"— 



or magnetii 



,39.40 



has permitted the realization of one- 



dimensional systems of interacting particles. In paral- 
lel, mixtures of bosonic and fermionic atoms^i^ het- 
eronuclear mixtures of fermionic atoms^i^ as well as 
pseudospin-1/2 fermionic atoms^~ have been trapped 
and cooled. Recently, partially polarized pseudospin-1/2 
fermionic atoms have been trapped in an array of ID 
tubes^. Finally, an analog of of photoemission spec- 
troscopy for cold atomic gases has been developed^li^ 
These developments might permit the future measure- 
ment of spectral functions of one-dimensional trapped 
atomic gases and a comparison with theoretical predic- 
tions. 

From the theoretical point of view, spectral functions 
A{q,uj) for the TL liquid at T = have been consid- 
ered in Refs. 14911501 . Power-law singularities with ex- 
ponents depending on the TL-liquid parameters have 
been predicted at oj = ±Mi/|q| {v = p,<j). At T > 
0, the spectral function, A{q,ijj), is strongly affected 
by the thermal fluctuation, which reduces the effect of 
interaction. Actually, double peaks due to the spin 
charge separation move to a single peak with increas- 
ing temperature.'^^ The spectral function of the Hubbard 
model in one dimension have been investigated in the 
limit of [/ -> -i-oo^i^ with the help of the Ogata-Shiba 
wavefunction.^"* More recently, combining a bosonization 



approach for the charge degree of freedom and an ex- 
act treatment of the spin degrees of freedom has per- 
mitted to obtain an analytic expression of the spectral 
function in that limit™. For finite U, the Hubbard model 
has been considered using exact diagonalizations^^, quan- 
tum Monte Carlo,"^^ and dynamical density-matrix- 
renormalization-group methods .^^'^^ The spectral func- 
tions of a magnetized TL liquid have been studied 
analyticallySii^ as well as numericallyi^ 

In the present paper, we calculate the fermion spectral 
function in a two-component TL liquid for zero temper- 
ature. In Sec. ini we recall the bosonization treatment 
of the two-band of spinless fermion modeljM and derive 
the expression of the real-space fermion Green's func- 
tion at zero temperature. The bosonized Hamiltonian 
of the two-band spinless fermion models also describes 
spin-1/2 fermions in a magnetic field^^— and mixtures 
of spinless fermions with bosons or fermion o^^'^^ . In Sec. 
llVi the spectral functions are obtained. In the case of 
an SU(2) invariant model, corresponding to a TL liquid 
of spin-1/2 fermions in the absence of magnetic field, the 
fermion spectral function can be expressed in terms of the 
Gauss hypergeometric function as shown in Ref.uOlln the 
non SU(2) invariant case, which might be achieved exper- 
imentally with polarized spin-1/2 neutral fermions,— or 
quantum wires under a magnetic fleld^i ^jj^g formion spec- 
tral functions can be expressed in terms of Appell hyper- 
geometric functions. Our approach recovers the previous 
results,—'^ but also allows to describe the behavior of 
the spectral function away from the points uj = ±u^q. 
We discuss the applications of our results to the Hub- 
bard model under a magnetic field in Sec. IIV B 31 



II. BOSONIZATION 

We consider a general two-band model of interacting 
spinless fermions defined by: 



H 



a = l,2 



dxVa{^lj]i^a^x1pRM " ^l ^a^x^l^ L ,a) 



+ E / ^^Sapl +9 dxpip2, (2.1) 

where ^pR^a and ipL,a respectively annihilate one right 
moving and left moving fermion in band a, pa — 
'^R,a'^R:a + i'L,a'^L,a, Va is thc vclocity of fermious in 
band a € {1,2} and ga the strength of intcrband in- 
teraction. The model includes only the forward scat- 
tering interaction and not the backward scattering ~ 
"^R i^L iV^L 2^R 2- This assumption is justified when the 
two bands have different Fermi wavevectors or when the 
backward scattering interactions are irrelevant. 

The model (|2.ip can be bosonized perturbatively, 
yielding: 



H = 



E 

a=l,2 



dx 
2^ 






-^ I dxdx4>idx4>2, 



(2.2) 



with [(/>a(x), Hf,(x')] = iSa,bS{x — x') and a,b £ {1, 2}. In 
Eq. (pr^ . we have: 



"a = «a Va 



'29a 



Ka ^ [I 



'29a 



-1/2 



(2.3a) 
(2.3b) 



The intraband interaction is included in Ki and K2 ■ The 
case Ka < 1 {Ka > 1) corresponds to the repulsive (at- 
tractive) interaction. 

More general models (for instance lattice models) can 
also be considered with bosonization. A non-perturbative 
formulation leads to a Hamiltonian^ 

/J 
^ [TT^MabU^nk + N,bdx<Padxh] , (2.4) 

a,o 

where the matrices M and N are real symmetric and 
are defined in terms of the variations of the ground state 
energy Eqs of a finite system of size L from (respectively) 
change of boundary conditions ipa{L) = e''^"'0a(O) and 
change of particle densities pa — Na/L: 



Mab = ttL 

Nab 



difadipb ' 

1 O^Egs 



ttL dpadpb 



(2.5) 
(2.6) 



The spectrum of the general bosonized Hamilto- 
nian (j2.4p is obtained by a linear transformation of the 
fields Ha and (j)a- 



Hfa - J^Pbpflf}. 



,Qa 



(2.7) 



(2.8) 



where P*Q = 1 in order to preserve the canonical commu- 
tation relations. ^^ The matrices P and Q are calculated 
explicitly by applying a succession of linear transforma- 
tions. First, the matrix M is diagonalized by a rotation 
Ri (Ai = *RiMRi) which transforms the matrix TV into 
Ni — *RiNRi. Hereafter we denote <f) — *{(f>i,4'2) and 
n ~ *(Hi,n2). By the transformation 11 — i?ini and 
(p = Ri4>i, the Hamiltonian is thus transformed into: 

r dr 

H = J ^ [TT^'niAiHi + *(9,0i)Ari(a,</.i)] . (2.9) 

— 1/2 

Using a second transformation IIi = A^ II2 and (pi — 
Aj^' 02, the Hamiltonian becomes: 

dx 



^='2;^ 



^2*n2n2 + *(a,02)Ai/'iViAi/'(9,02) 



1/2/ 



(2.10) 



As the matrix Aj^' NiA^' is symmetric, it can be diag- 

onalized by a rotation i?2 i.e., A^^ NiA^ — i?2A2*i?2- 
Writing 112 = -^2113 and 02 = R2<p3, we find that the 
Hamihonian can be diagonahzed as 

i? = y ^ [vr^ 'HsHa + *(a,03)A2(a.03)] , (2.11) 

in which the modes are decoupled. Finally, we can rescale 
the fields n3 = (Aa)^/''!! and 03 = (A2)-i/''0 to write: 



H= 1^ 

2tt 



n^ *fl(A2)i/2fl + *(a,0)(A2)i/2(a,0)' 



(2.12) 
In this last equation, the elements on the diagonal of 



(A2)^/^ are the velocities up of the decoupled modes of 
the Hamiltonian (|2.4p . The stability of the multicompo- 
nent TL liquid state requires that the velocities in (j2.1ip 
are real, i.e., that the matrix MN has only positive eigen- 
values. The transformations can be written explicitly as: 



p = i^lA-l/'i?2(A2)l/^ 

Q = i^lA}/'i?2(A2)-l/^ 



(2.13) 
(2.14) 



and we have: *PAfP = (A2)^/2 and ^QNQ ^ 
(A2)^/^. This implies in particular that: ^PMNQ = A2 
i.e. Q~^MNQ = A2, and by taking the transpose, 
P^^NMP — A2. The excitation velocities u± where 
A2 = diag(u'^,w^) are obtained as 



4 



M12N12 ± 



{M11N12 + M12N22) {M21N11 + M22N21), 




(2.15) 



r 



where u\ = NuMn and 1*2 = N22M22- 

The diagonalization method allows us to derive also 
expressions for the Green's functions of the chiral fields. 



•JRm 



i>L,a 



(2.16) 



where 9a ^ t^ J dx'IIa{x'). The field operators are ex- 
pressed as 



V'fl,a(a;,i) 
V'L,a(a;,i) 



1 



'27ra 

1 
e 

'27ra 



'i(t>n^a.{x,t) 



.i<t'L.a(x,t) 



(2.17a) 
(2.17b) 



where a is the short-distance cutoff. The Hamiltonian 
(|2.12[) can be reexpressed in terms of non-interacting chi- 
ral fields: 



H 



dx r 

47r 



.1/2. 



\d.,ct>R)A'2''{d^CJ>R 

with the transformation: 

(t>L = I \{Q - P)4>R. + {Q + P)4>L 



\dM^y\dx4>L) 



(2.18) 

(2.19) 
(2.20) 



where 0_r = \4'b.,iAr,2) and 0l = \4>L,i,(t>L,2)- 

For the model (|2.ip , the off-diagonal term of N is given 
by A'^12 — g /t^ and there is no interband current-current 
interaction, i.e., M12 = 0. In this case, the explicit forms 
of the matrices P and Q for the model (|2.ip can be ex- 
pressed in a compact form. Since M12 = 0, the matrix Ri 



becomes unit matrix, and Ai ~ M — diag(Mi_R'i,'U2^2) 
The matrix A2 is given by A2 — diag(u^, u^) where 



2 _ul + u^ 



± 



U1K1U2K2 



2 V V 2 y \7r/ 

(2.21) 
We note that the velocities u± depend on g^, i.e., do not 
depend on the sign of g. From Eq. (|2.13p . we obtain 



P 



Q 




where tana = 2{g / tt)\/ U1K1U2K2 / {u\ ~ u\). From Eq. 
(|2.2ip . the stability condition is given by 



U1U2 > 



{^yK.K2. 



(2.24) 



In the following analysis, we restrict ourselves to the case 
of ui > U2, and Eq. (I2.24p gives the lower condition for 
M2, i.e., W2 > U2c = {g/Tr)'^KiK2/ui. 



III. PHASE DIAGRAM 

We derive a phase diagram for the model (|2.ip . by 
examining asymptotic behavior of correlation functions: 



{Oa{x)Oa{0))^x-'^^, 



(3.1) 



where the Oa's represent the order parameters and the 
r^yi's the correspondmg exponents. We restrict ourselves 
to the simplest case gi — 52 in Eq. ()2.1|) . and assume 
1'2/wi < 1 without loss of generality. In the TL-liquid 
state, the state with the smallest exponent represents the 
(quasi) long-range-ordered state and thus we can deter- 
mine the phase diagram. As possible order parameters, 
we can consider the intraband density wave (DW) and 
pairing of superconducting state (SC), which are given 
by Oowa « exp{i2(j>a) and Osca oc exp{i29a). Their 
respective exponents rjA are given by 



jyDWi 



'7DW2 



?7sci 



?7SC2 



uiKi 
U2K2 

U+U- 
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1 
U2K2 
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ul 
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ul 
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uV 
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(3.2a) 
(3.2b) 
(3.2c) 
(3.2d) 



In addition, we can also consider the order param- 
eters for the interband DW and interband SC states, 
given by OintorbandDW = ipR.i'^La oc e*"*«.i+*'^^-2 and 

Ointcrbandsc = ^R,i^L,2 « e-">«.i+'*^.= . The corrc- 
spending exponents are given by 



?yintcrband DW 



'7intcrband SC 



(gM 



u+ 

7?DW1 



U- 
'7DW2 



??DW1 r]DW2 ??SC1 ^SC2 

4 4 4 4 

U1K1U2K2 

U+U- 
?/SCl ??SC2 

4 4 



{9/^) 



1 



, (3.3a) 



, (3.3b) 



In Fig. [TJ these exponents are shown as a function of 
V2IV1, for repulsive interactions. For ^2/^1 = 1, the in- 
terband DW state becomes dominant \i g > 2ga [see 
Fig. [T] (a)] while the DW2 state becomes dominant if 
g < 2ga- For decreasing V2, the interband DW state is 
unfavorable and instead the DW2 state becomes domi- 
nant since the effect of the intraband interaction ga/v2 
is enhanced. On the other hand, if the intraband inter- 
actions are absent {ga = 0) [see Fig. [T] (b)], the DW2 
state is no longer enhanced for small V2 and instead the 
exponent 77302 decreases with decreasing V2 and the SC2 
state becomes most dominant state and finally the two- 
component TLL state becomes unstable {u^L < 0). In 
Fig. [21 the phase diagram on the plane of ^2/^1 and 
gal {t^vi) is shown with fixed g/ {t:vi) = 0.3. In the anal- 
ogy to the spinful electron model, the DW2, interband 
DW, and SC2 states corresponds to the conventional 
charge-density-wave, spin-density-wave, and triplet SC 
states, respectively. The region of the interband DW 
state becomes narrow with decreasing V2 and shrinks at 
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FIG. 1: (Color online) Exponents via for possible order pa- 
rameters as a function of V2/V1 for ga/(nvi) =0.1 (a) and 
ga/iiTvi) = 0.0 (b) with fixed g/{nvi) — 0.3. 



the point iv2/vi,ga/{nvi)) sa (0,0.04), where the SC2, 
DW2 and interband DW states have the same exponent. 
We note that for g = 0, we can recover the behavior 
that the DW2 state becomes dominant for ga > while 
the SC2 state becomes dominant for — 7rw2/2 < ga < 0, 
and the system becomes unstable if ga < —ttv2/2. We 
note that the expressions (|2.3p are valid only for weak 
interactions, and thus the precise determination of phase 
diagram for the region ga/v2 ^ 1 in Fig. [2] is beyond the 
present approach. From the qualitative considerations, 
the following modifications to the phase diagram can be 
expected. Since the parameter K2 would take a nonzero 
value for ga/v2 — >■ 00, the condition (|2.24p cannot be 
fulfilled for small W2/W1, and then the unstable region 
always appears in the small limit of W2/f 1 for all values 
of ga. Furthermore, by noting that the unstable region 
is adjacent to the SC2 state, it is expected that the SC2 
state is obtained even for large ga region and is located 
in between the DW2 state and the unstable region. 



IV. SPECTRAL FUNCTION 



The spectral function is obtained from 



1 



Aa{k,Lj) = Im / dxdte-'^^''-'^'^Ga{x,t), (4.1) 



where Ga{x,t) is the retarded fermion Green's function 
Gaix,t) = -i6'(i)({V'a(a;,i),V'l(0,0)}),. Introducing the 
right moving and left moving components, 'ipa{x^t) = 



0.2 



DW2 




FIG. 2: (Color online) Phase diagram on a plane of v-^jvx and 
ga/iiTVi) with fixed g/ijivi) — 0.3. In the "unstable" region, 
U2 becomes imaginary. 



J 



e*''^'"^Vfl:,a(2:,i) + e '■'"'•'''^tpL,a{x,t), the chiral Green's 
functions {v = R, L) are given by 

where F,y_a{x,t) — {^l;^_a{x,t)'^pl^^{0,0)). Then the spec- 
tral function can be decomposed into the four contribu- 
tions: 

A^^a{q,^) = /i.,Q(g,a;) +/t,,Q(-g, -w), (4.3) 

where 



/^,a(g,w) = 



27r 



rfte*"* 



dxe-"'''F,,a{x,t). (4.4) 



The direct calculation of the Green's functions for the 
phase variables yields 



((61a (a;, t) - 4>a{x, t)){eaiO, 0) - MO, O)))conn. = J^ "'^'^^ ^^ 



;3=± 



{{ea{x,t) + Mx,t)){0aiO,O) + MO,0)))conn. = ^ </3 In 

where the exponents i^s are given by 



a + i{upt — x) 

a 
a + i{upt — x) 



'^Lfi In 



'^a,B In 



a + i{upt + x) 
a 



a + i{upt + x) 



Va..li = -riPaf! + Qapf: v'a,fi = jC^a/J " Q a/if ■ 



, (4.5) 
, (4.6) 

(4.7) 



Thus the one-particle Green's functions are expressed as 



(V'fl,a(x,i)^lj,a(0,0)) = ;^n 



2Tra 



(V'L,a(x,i)V'L(0,0)) = ;^n 



27ra 



a + i{uijt — x) 



a 
a + i{upt — x) 



a -t- i{uijt + x) 



a 
a + i{upt + x) 



(4.8a) 
(4.8b) 



Substituting the expression from Eq. (|4.8p into Eq. (|4.4I) . we obtain the integral form for the spectral function. We 
note that we have the identity: X^fll'^a,^ ~ ^'a b) ~ 1- '^^^ spectral function is obtained from the integral: 



IrAi,^) 



(2^)2 



dx dt- 



pi(ujt — qx) 



[a + i{u+t - x)\''-.+ \a + i{u^t - x)]"'^.- [a + i{u+t + x)Y''-+ [a + i{u^t + x)y---- 

I 



, (4.9) 



where Va = {va,+ + Va,-+v'a + + v'g^.~)- A similar expres- 
sion for left moving fermions with i'a,± and v'^^ ^ inter- 
changed. 



A. SU(2) symmetric model 



Let us first consider the case with SU{2) symmetry and 
repulsive interactions, and show how the expressions in 
terms of Gauss hypergeometric functions are recoveredi2£ 



The bosonized Hamiltonian reads H = Hn + H„ where 



Hp = 



H„ = — 



dx 
2^ 
dx 
2^ 
2gi± 
(27ra) 



Kp J 



dxcos 



(4.10) 



(4.11) 



with Up > u„, Kp < 1 and K^ > 1 for repulsive interac- 
tions. Under the renormahzation group, K^,gi± flow to 
fixed point values K* = l,gl± = 0. The resulting fixed 
point Hamiltonian is therefore in the diagonalized form 
of Eq. (|2.12p . We will approximate the spectral function 
of this model by replacing the exact Green's function by 
its fixed point value. This amounts to neglect logarith- 
mic corrections. Within this approximation, the spectral 
function of right-moving fermions of spin s —"[, I is given 
by AR^s{q,i^) = lR{q,^) + lR{-q,-Lo) with 



iniq,^) 



v27p 



(2^)2 



dxdt e 



i{ujt~qx) 



+ i[upt - x)]'>'p+i/2 [a -t- i{u„t - a;)]i/2 



+ i(upf + x)]>' 



(4.12) 



where 7p — [Kp + K ^ — 2)/8.— The correspondence 
between Eqs. ()4.9p and (|4.12p is given by m+ = Up, u_ — 
Ua, i'a,+ = (7p + l/2), Ua- = 1/2, i^a.+ =Jp, and iy'^_ = 
0. To calculate the integral in Eq. ()4.12p . we use the 
Feynman representation^^: 



1 



r(z^i + V2) 



,,1^1-1 



dw- 



{1-w) 



1^2 — 1 



(4.13) 
which is valid for i^i, i^2 > 0. Then, with the help of (jA4l) . 
we obtain: 



Iniq,^) 



nip + 1) 



^7p-1/2(1_^)-1/2q,7p 



dxdte^'^"*''^'-'^ / dv 

4^2r(7^ + l/2)r(l/2) 7 Jo {a + t[vup + {1 - v)ua]t - ix}''p+^ 



a + i{upt + x) 



ip 



After the space-time integration, we have: 

iR{q,^) 



r(7p + i/2)r(7p)r(i/2) 

1^ - [vup + (1 - v)u„)q\'^''^'^ 



lu + Upqp" dw>-l/2(l-i;)-l/2 



Q{uj + Upq)Q{uj - [vUp + (1 - v)ua]q), 



r 



(4.14) 



(4.15) 



where 9(x) is the Heaviside step function, and r{z) is the 
Euler Gamma function. We will consider the case oi q > 
0. From Eq. (|4.15p . we find Iniq, uj) = ioi uj < u^q and 



lR{-q,- 



for UJ > —Upq. Therefore, Aji^iqi^) — 



when —Upq < lu < u„q, while Aji_s{q,uj) — Iji{q,uj) for 
UJ > Ua-q and Afl^s(a,a;) = lR{—q,—Uj) for uj < —Upq. 
For the calculation of /_r((7, uj) when uj > Ua-q we have to 
separate the two cases, u^q < uj < Upq and uj > Upq. 



1. U} > UpQ 



For UJ > Upq, the two Theta functions in Eq. (|4.15l) can 
be replaced by one. With the change of variable 



w 



2upV 



Up + u„ + {up -u„)v' 



the integral Eq. (|4.15p becomes: 

Ad ^(q.uj) = — ; — -, 

^'^^''' ' r(7p)r(7p + i) 

{uj + Upq)^" {uj - u^q)''"^'^ 

^ (2Mp)^P + l/2(y^ + U^)7p-l/2 

X 2^1 ( l-7p,7p+ o^^p + lr 



Ua UJ -\- Upq 



2Ur, 



Uaq 



(4.17) 



For w — !► oo, I{q,uj) ^ uj'^"''' ^. The expression ()4.17p 
agrees with Eq. (19.27) of Ref. [t^ with the notation 



6* = 27p. When uj 



+ 0, the argument of the hy- 



pergeometric function in Eq. ()4.17p becomes equal to one 
and for 7p < 1/2, the hypergeometric function is diver- 
gent. Using Eq. (jA6p . we rewrite (I4.17P as: 



(4.16) AR,,{q,uj) 



(a/2up)2> 

r(7p)r(7p + 1) 



(a; + Mpg)>(a;-Upg)'^''-i/2 



{uj - u^qy/'^ 



X 2F1 (27p,-;7p 



Up — Ua UJ + UpQ 

2up u! — u^qj 

(4.18) 



in which the hypergeometric function remains finite 
as w — 7- Upq + 0. We then find that Ar_s{<1,^) ^ 

(2«p)^p(«p-u„)i/2r(l/2+7p)7r y^ ^Pl) as UJ r^ Upq. 

The power-law divergence was previously obtained by 
Voit^- and by Meden and Schocnhammer'*^ by analyz- 
ing the divergence of the integral (|4.12p in the vicinity 
of w ~ Upq. The expression in terms of hypergeometric 
functions also provides the prefactors. 
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Ua-q < IJJ < UpQ 



When Ua-q < uj < Upq, the integration over v in ()4.15p 
is limited to the range < v < {uj — Uaq)/{upq — u^q). 
With the change of variable: 



w 



Upq 



(up - u„)v 



UJ - U^q [Up + Ua) + {Up - Ua)v ' 

the integral (|4.15p reduces to: 

a27p 



(4.19) 



lfl,,s 



(9,w) 



r(l/2)r(27p -M/2) 



{uj + Wpq)-l/2(^ _ y^^)27p-l/2 
(Up-fU,)7p-l/2(up-U,)7P + l/2 



X2i^l ( -,lp + -]2-fp 



Uaq 



2 Up — Ua uj + Upq 



(4.20) 



For u! — ?> u^q + 0, the expression (|4.20p has a power law 
divergencej_~_(ci;— Mo-q)^'''''^^/^ for7p < 1/4, in agreement 
with Refs l49li50l . When uj — >■ Upq — 0, the argument of the 
hypergeometric function becomes equal to one, leading 
to a power law divergence. Using again Eq. (jA6p , we can 
rewrite the Eq. (|4.20p as: 



^i?„s(g,w) 



v27p 



r(i/2)r(27p + i/2) 

(uj ~ Uaqf^''^'^''^{upq - uj)'^"-'^/'^ 
{uj + Upq)''p{up-Ua)^^p 

2u„ UJ - 



X2i^i 27p,7p;27p4- 



u^g 



2' Up — Ucr W + UpQ 



and recover the divergenci 

a2^Pg"'p-V2r(l/2-7p) 



.49.50 



C^)7p-1/2 



(4.21) 



for 



(2«p)"'P(«p-«„)i/2r(l/2+7p)7r("P? 

UJ — >■ Up(7 — 0. Although the exponent is the same 
on both sides of a; = Upq, the peak is asymmetric, the 
ratio of amplitudes being sin(7r7p). 



FIG. 3: (Color online) The spectral functions Aii,s{q = 
0.01, a;) with SU(2) symmetry for several choices of 7p (= 0.2, 
0.4, 0.6, 0.8 from top to bottom) with fixed Up = 2, Uo- = 1, 
q — 0.01. The dots represent the value at a; = Upq given by 
Eq. gSll). 



3. UJ < —Upq 

In the case uj < —Upq, we need lR,cr{—q., —uj). In the 
integral (|4.15p . both Theta functions are then equal to 
one, and with the change of variable: 



2upV 



{Up + Ua) + {Up ~ Ua)v ' 



(4.22) 



we find: 

AR.a{q,i^) 



a^^p 



F(7p)F(7p + 1) 
\uj - Uaq\'^p~'^\u} + Upql'^p 

(Up + U,)7p-l/2(2up)7p + l/2 



X 2F1 [ i~ipnp + ^;ip 



Ua UJ + Upq 



2Ur. 



UJ — Uaq 



(4.23) 



For UJ — > 
ishes with 



—Upq — 0, the expression (|4.23p van- 
a cusp singularity Aji^a{—(1,—^) ^ 

'''' , in agreement 



l^TpoTp-l 



(2np)^Pr(7p)rf7„-l-.l)[2Mp(«p+«„)]i/2 

with Refs. 49.50;. 



\UJ- 



-Upq\ 



Typical behavior of the spectral functions A/j ,j(g.cj) 
is shown in Fig. [S] There is no weight at —Upq < uj < 
Uaq. The overall profile reproduces the previous results 
given in Refs. ligllsol . We note that for 7p > 1/2, the 
divergence at a; = Mpg disappears and is replaced by the 
cusp structure. From Eqs. (|4.17p and (I4.20p . its peak 
value is given by 



AR,s{q,Upq) 



a^>r(7p - ^) (Au)7p-lg27p-l 

r(7p)r(27p)r(i)(2up)i/2(2u)7P-i/2 



(4.24) 



where u = {up + Ucr)/'2. and Au — [up — Ua)- In Fig. |3l and the asymptotic behavior for i < 7p < 1 is given by 
the peak positions are represented by the dots. Thus the 
asymptotic behavior at a; w Upq for 7p < ^ is given by 



AR,s{q,^) 



C{uj — Upq)^p 2 sin7r7p (cj > Upq) 



C{upq~uj)''^- 



iu) < Upq) 

(4.25) 



AR,s{q,(^) 



AR^siq, Upq) - C"(a; - Upq)'^" 2 sin7r7p (u > Upq) 

AFi,s{q, Upq) - C'(upq - uj)''''^^ {u < Upq) 
I 



(4.26) 



where C and C" are positive numerical constants depend- 
ing on q. Because of our simphfied treatment of the cut- 
off, the sum rule J dijjA{q, lj) — 1 cannot be satisfied. 
In a more rigorous treatment, the short-range cutoff a 
used in the construction of the creation and annihilation 
operators^i and the momentum cutoff for the interactions 
must be treated independently!^ However, in our paper, 
we are only concerned with the asymptotic behavior of 
the spectral functions for momenta that deviate from the 
Fermi momenta by an amount which is much less that 
the momentum cutoff. In such a case, the corrections 



J 



resulting from having two-distinct cutoffs can be safely 
ignored. 



B. General two-band model 

Next we consider the case of general two-band model, 
and evaluate the integral (|4.9p . By using twice the Feyn- 
man identity (|4.13|) and Eq. ()A4[) . we can reexpress Eq. 
(gH) as 



lR,a{q^^) 



^,i^a-l 



ri,ya,+)niya.-)T{iy',^+)riiy',^_) Jo jo 

x[uj - u{wi)q]<- + +<-~^[uj + u{w2)qV-^+''' 

I 



dwi I d'W2u/['^'^ (1 — lui )"'"'" ^w^ 



..+ -i(l_^,)<--i 



_^i Q[u) — u{wi)q\Q[ui + u{w2)q\ 
[u{wi) + u{w2)Y--'^ 



(4.27) 



where u(w) = wu+ -I- (1 — w)u-. From Eq. (|4.27p . we 
obtain A^^aiq > 0, |a;| < U-q) = 0. Contrarily to the 
SU(2) symmetric case of Sec. IIVAI the spectral function 
does not vanish anymore when —u^q < uj < —u^q. 



1. Power-law singularities 

The power-law singularitie s '^^1^° can be recovered from 
Eq. (|4.27l) in a simple manner. In the analysis of the 
power-law singularities, we restrict ourselves to the case 
g > 0. Indeed, if we consider the case of w — )■ u+q + 0, 
the dominant contribution to the integral comes from the 
integration over wi in the vicinity of wi = 1. The integral 
to consider is then: 

lR,a{q,^) oc / tiwi(l - wi)''"'-^^ 
^0 



X [(w - u+q) + Aug(l - wi)Y-^--+'''^.+-^ 



\uj-u+q\ 



-+y' .+y' +-1 



|A^g|-^" 



(4.28) 



provided that Vg,.- + '^a - + '^a + < Ij where Au = 
(m+ — u_). Similarly, when u^q < iv < u+q, the 
wi integration also determines the power law singular- 
ities. This time, the integration over wi is restricted 
to < wi < {uj — U-q)/{Auq). Changing variables to 
Hji — wi{uj — u^q)/{Auq), we have to consider the inte- 
gral: 



(a; - u^q)'' 



(Am(j)''".+ -i 
X {1 - wi)<-+'''^-+-^ 



dwiWi" 



UJ — U-q 
Auq 



(4.29) 



When uj -^ u^q + 0, the integral goes to a con- 
stant, and we have the power law; Ijiail,'-^) ^ {io — 
M_g)''°++''-++''»-"^ When uj -^ u+q - 0, the in- 
tegral has a power law singularity, and the behavior 
lR^a{q,^) ^ (-W + u+q)''°--+''».++''"-"^ is obtained. 
Such behavior was previously obtained in Refs. 49il50l. 
By a similar method, one can also obtain the power law 
singularities at w = —u±q of lR^a{—q,—to). This time, 
the origin of the singularities is the integration over 'W2 ■ 
By summarizing, the asymptotic behavior of the spectral 
function are given by 



AR^aiq,^) OC 



where 



f 


OJ — u+q^°-+ (for w ^- +M_|_q ± 0) 




{u} — u^q)^"-- (for w ^- +M_g + 0) 




{uj + u^q)'^'^-- (for cj ^ - 


-u_g-0) 




C +\uj + u+qf'^-+ (for UJ ^ - 


-u+q±Q) 
(4.30) 


Pa 


+ = >^a- + v'^.+ + v'^- - 1, 


(4.31a) 


Pa 


- = Va^+ + v'a^+ + v'^,_ - 1, 


(4.31b) 


P'a 


- = Va,+ + Va- + v'a^+ - 1, 


(4.31c) 


P'a 


+ = Va^+ + Va- + l^'a.- - 1, 


(4.31d) 



As Ua,+ + Va- = 1 + i'a+ + t^'a -: there is no divergence 
but only a cusp in the vicinity of a; = —u±q. 

From Eq. (14.71) . the exponents which determine 
the Green's function for the right moving particle in 
Eqs. (|4.9p are given by 



'^i.± 



i,± 



J^2,± 



^2,± 



1 fKiUi 

u± 
Kiui 

u± 

K2U2 

u± 

K2U2 

u± 



KlUi 

u± 
KlUi 

u± 
K2U2 

u± 
K2U2 



+ 2 
-2 
+ 2 
-2 



1± 



1± 



It 



It 



ul^ 


^ul 


ul- 


-ul 


u\- 


-u^_ 


u\- 


-ul 


u\- 


-v?_ 


u\- 


-ul 


<- 


-W2_ 


a. (14.321) 



(4.32) 



Figure |4] shows the exponents Va^p, [Eq. (I4.32p ] as a 
function of difference in velocities. In order to distin- 
guish the effects of the intraband coupling g^ and the 
interband coupling g, we plot the U2/U1 dependence of 
the exponents with fixed Ka (i.e., fixed Qa/va)- In this 
case, the two-component TL liquid is always unstable 
for small U2 , in contrast to the situation in Fig. [51 For 
Ml = U2, the behavior is similar to that of the TL liq- 
uid with only the intraband forward scattering where 
(i^a.+,!^a,-, <+,<_) « (0.514,0.501,0.014,0.001), sug- 
gesting the weak effect of the intraband interactions. 
With decreasing U2, the exponent vi^^ (j^2,+) decreases 
suggesting that the integral at a; « u-t [x ~ u^t) for Ir^i 
{Ir,2) in Eq. (j4.9p becomes less singular. Although the 
effect of intraband interaction on v'l ^ and 1^2 ± is small 
for U2/U1 ~ 1, the exponents v'l „ and 1^2 - are enhanced 
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FIG. 4: The U2/u\ dependence of the exponents /^i,±, u[ ^ 
(top) and 1^2, ±, 2^2, ± (bottom) for K\ — K2 = 0.8, and 
g/{TTUi) = 0.3. The two-component TL liquid is unstable 
for U2 < U2c ~ 0.0576mi. 




0.5 
U2/U1 



FIG. 5: The U2/u\ dependence of the exponents /?i,±, Pi,± 
(top) and /32,±, /?2,± (bottom) for K\ = K2 = 0.8, and 
g/{-KUi) = 0.3. 



for small U2 just above U2c- In Fig. [Sj the correspond- 
ing Pi,±, fi2,±, P'l ±, and j32 j. are shown. The case of 
K = 1 (not shown) is similar to Fig. [5] Note that the 
exponent v'l j_ becomes extremely small but nonzero for 
Ki = K2 — 1 and small .g/vr. Its asymptotic behavior is 
given by iy[^^ ~ [g/{TTUi)]*{u2/ui)'^ /16 for U2 < ui. 

As will be shown in the next section, the spectral func- 
tion has two divergences at w = u+q and uj = u^q when 
fi2,+ < and /3i,_ < 0, respectively. However each diver- 
gence is replaced by a cusp when (32,+ > and /3i,_ > 0, 
respectively. Figure [H] shows such a critical value as a 
function of g/{TTUi) with fixed /32,+ = (solid lines) and 
/3i__ = (dotted lines) for two cases where the intra- 
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FIG. 6: The values for U2/U1 satisfying /32,+ = (solid line) 
and Pi- = (dotted line), with fixed K{^ Ki = K2) = 1 
and 0.8. 



band interaction is absent {Ki = K2 = 1) and is present 
{Ki — K2 — 0.8). For Ka = 1, the critical value of U2 in- 
creases as a function of g while that takes a minimum for 



K = 0.8. The minimum for if ^ 1 suggests a competi- 
tion between the intraband interaction and the interband 
interaction. With increasing K tol, the value of g/{nui) 
at which U2 takes a minimum decreases, and the line in 
the limit oi K -^ 1 coincides with that of if = 1 except 
for g/iT — 0, where M2/U1 = 1 at g/n = is always ob- 
tained for if y^ 0. Thus it is found that double peaks of 
either 1 or 2 particle is replaced by a single peak in the 
region enclosed by the solid line and the dotted line in 
Fig. El 



2. Representation of spectral function as integrals of 
hypergeometric functions 



It is possible to rewrite the double integral (|4.27p as 
a single integral containing the Gauss hypergeometric 
function^ 2i^i(a, /3; 7; z) (see Appendix [X)) . We again 
restrict ourselves to the case of g > 0. The spec- 
tral function for q < can easily obtained by noting 
A^^aiq, oj) = A^,a(-9, -oj) [see Eq. (gSl)]. For |a;| > u+q, 
we obtain 



J 



{qay-- 



\VcL,-+y'„ ++'^'n _-l 



(w - 'a+g)'^--+'^-++''— -^(w + u+qY-^-+ 



l^a, + +'^a.-—l 



AnA<l,-^)Uu,, r(.„,+ + .,,_)r«+)r«_) (,.-._g)^.-(2u+g)--++-.-+<.-i(2ug)<,- 



X / dt {l-t)<'+-^{t)<--' 



^fl,a(g,^)L<-« 



+9 



^ ^ Au{uj-u+q) ^ 
2u{uj + u+q) 

X 2i^l {Va - 1, l^a.-; Va^+ + Va.^\ Zl) , (4.33a) 

jqc^f--' {H - u+qr^^^+^^-+<--Wu:\ + u+q)<++'^'^--^ 

n<^+ + <,-)n^a,+)n^a,-) {\uj\ ~ u^q)<- {2u+qy-'++<++<--\2uqy^.~ 



X / dt {i-ty-'+'^{tY-'—^ 



^ _^ Au(H_u+qO ^ 



2m(|w| + u+q) 
X 2i^i (i^a - 1, v'^^--, v'^^+ + v'^^-\ zi) . 

Here u = {u+ + u_)/2 and Am = {u+ — uJ). The parameter zi is given by 



zi 



For u-q < \u)\ < u+q, we have 

^i?,a(<Z,'^)U_g<Lj<M+(j ^ 



Au(|aj| + u+q) 
2M_|_(|a;| — U-q) 



[qaf 



^ j ^u{\uj\ -u+q) ^ 
2u(|a;| + u+q) 



(4.33b) 



(4.34) 



AR,a{q,(^)\-u+q< 



LoK—U- q 



r(^a.+ + y[.,+ + u'aZT{VaZB{v',^+,ul_) 

(w - u_q)''"++''»> ++<.--l(-w + u+q)'''^'-+<-+< + -^ 
(w -I- u+q)''-+ {bJ + u^q)<.- (AMg)P"-i 

X f dtil- t)'''-^-\t)<-~' 2F1 {Da " 1, iya.+ + 1^'a,-; ^a,+ + y'a.+ + y'a,--. Z2) , (4.35a) 

Jq 

(qa)''--^ 

r(l/a.+ + l^a- + y'a,+ )T{v'^_)B{Va^+,ya,A 
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{\uj\-u_q) 



Va, + +Va.-+'^'a.+ ^^( _ 



+ '{-\uj\+u+q) 



l^a.4--\'^a.— ~\~^n - ~^ 



{\uj\ + u+qY-:+ {\uj\ + u^qY--- (Auq)^--^ 

'0 



<+;z2), (4.35b) 



r 



where B{p,q) is the beta function B{p,q) 
T(p)T{q)/T{p + q). The parameter Z2 is given by 



Z2 



2uj^{\uj\ — u^q) 
Au(|a;| + u+q) 



^ ^ l\u{-\uj\ + u+q) ^ 
2u+{\(jj\ + u-q) 



(4.36) 



We find AR,a{q,^)\-u^q<uj<U-q = 0. 

We note that, if v'^_ = 0, we have Aji^a — for 
—u^q < u < —u^q < 0. In Eq. (j4.33p . the power-law 



J 



singularities are already factored out and the integrals 
are regular. This fact allows the numerical evaluation of 
the formulae (|I351) and (|T^ . 



The integrals in Eq. ()4.33p can also be reduced to 
the Appell hypergeometric functions^SiIl. By using Eq. 
(jAlip . the spectral functions for the range U-q < |w| < 
u+q are expressed as 



AR,a{q,'^% 



(a/Au) 



i>a.-l 



{H-u.qY 



■+-'{-\Lu\+u+qy 



(n/AiiV 
AR^a{q,l^)\-u^q<uj<-u_q 



{\io\+u+qr^^+i\co\+u^qr^-- 
I 2u+{\ljj\ — U-q) 2u(|a;| — U-q) 



{a/AuY 



, (4 37) 
Au(|a;| +u+g) ' Au(|a;| +u-g)y ' 

{\uj\ - u_9)'^-++'^»-+''-+-i(-|w| + u+qY--++'''^^-+<--^ 



r(t^a,+ + ^a.- + <+)r«_) 

■X. Fiiva- YVa,+ ,Va-;Va,+ + >^a- + '^a,+ : 



i\Lj\+u+qY--+{\Lj\+U-qY-.- 
I 2u+{\uj\ — u^q) 2u{\uj\ — u^q) 
Au(|a;| + u+q) ' Au(|a;| + u^q) 



, (4.38) 



where -^1(0; /?, /?'; 7; a;, y) is the first Appell hypergeometric function of two variables. Similarly, by using Eq. (jA12p . 
equation (|4.35p can also be reduced to 



^fl.a(9,w)| 



{a/2u+Y 



iH -u+q) 



'^a.-+'^'a, + +'^'a. 



Wuj\+u+q) 



l'a.+ +l'a 



"^"^' r(z.,,+ + z.,.„)r(^^^+ + <_) {\oj\~-U-q)-'^-{\Lo\+U-q)<- 

Au{\u]\+u+q) Au{\u]\ — u+q) 



XF2 {va-Y Va,^, v'a-, Va,+ +l^a-, v'^^+ + v'^-] 



2u+{\ijj\ — u^q) ' 2u+(|w| + u^q) 



, (4.39) 



AR^a{q,uj)\ 



u:<—u+q 



(a/2M_| 



APa-1 



{\uj\ -U+q) 



l'a..-\-+Va.-+l''a,. 



Wuj\+u+q) 



'^a.-+V„ ,+V„ 



r«+ + <-)r(^a,+ + ^a,-) 



{\uj\-U-q)'''^:-[\Lo\+U-q)-^ 



„ / / / / Au(\u)\ + u+q) Au(\uj\ — u+q) , , 

XF2 [Va-Y y'a,-. ya^-\ <.+ + ^1,-, i^a.+ + ^a,-; TT^^ ^, . V, ^ ^ \ , 4.40 

' ' -r , 2u+{\u)\ ~ u^q) 2u+{\uj\ + u^q) ' 



r 



where F2(a; /3, /3';7,7'; x,?/) is the second Appell hyper- 
geometric functio n^^'^^ . These results are reminiscent of 
the one from Ref. UM, where the Fourier transform of the 
2kp component of the density-density correlation func- 
tion of a TL liquid with different spin and charge veloc- 
ities was shown to be expressible in terms of the Appell 
hypergeometric function of two variables. We note also 
that by setting u+ ~ Up, u_ — Ua, Va,+ = Yip + l/2)i 
Va,- = 1/2, v'^^+ = -ip, and v'^_ = 6 in Eqs. B371) . 
(1133), and (|OD)) . we can reproduce Eqs. (g^, (|il^ . 
and (|4.23p . respectively, showing that the Appell hyper- 



geometric function representation also covers the SU(2) 
invariant case. 

The special case of SU(2) symmetry can also be recov- 
ered from the integral representations (|4.33p and (|4.35p . 
However, to do this, we first need to consider the limit 
^^2 ^ in Eq. (|4T3)) . Rewriting (|4T3)) as: 



1 



_ r(t^i + V2) 

AYA-^- " r(j.i)r(^2 + i) 



i/i-i 



dw- 



i^2W'''~^(l -w) 

\AY\-w)^A2wY'^''' 



,(4.41) 
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FIG. 7: (Color online) Velocities as a function of the magnetic 
field for Ujt = 3 and p = 1/2. 



by making w — > (1 — w) and using V(y2 + 1) = t^2r(i^2) 
and changing integration variable to f = vf'^ , we find 
that: 



41^1 41^2 



r(;/i + V2) 

V{yx)Y{v2 
1 
dv- 



1 



(l_i,iM)^i-i 



[Ai(l-'yV''2) 4_^2wi/''2 



1'1+1'2 



(4.42) 

With this form, when 1/2 — > 0, v^/"^'^ — > OVu < 1, and the 
integral reduces to A^'^^ so that the Feynman identity 
remains applicable when 1^2 —^ 0. By applying the same 
transformations to Eq. (|4.33ap , we see that in the limit 
of i/^_ — )■ 0, the SU(2) case is recovered. The same re- 
sult also applies to Eq. (|4.35ap . In Eq. (|4.33b[) . when 
i'^_i_ — )• 0, the hypergeometric function in the integrand 
reduces to 1. The integral in Eq. (j4.33bp thus becomes 
a hypergeometric function 2^1 and the SU(2) invariant 
case is again recovered. Finally, in Eq. (j4.35bl) . a fac- 
tor r(f^ _) is present in the denominator of the fraction, 
but the integral remains finite in the limit ly'^ _ -^ 0. 
As a result, in the limit i^a + — ^ 0, the contribution of 
Eq. (|4.35bp to the spectral function vanishes, recovering 
fully the SU(2) invariant result. 



3. Application to the Hubbard model 

For the application of the present calculation to an ex- 
plicit model, here we consider the spin-i Hubbard model 
under the magnetic field: 

H = -t^(c] ,^Cj+i^^ + h.c.) -^^rij^^ 



-hj:( 



njA - nj,i) + UY^ rij^^rij.i, (4.43) 



where cr —"[, I refers to the spin degrees of freedom. We 
restrict ourselves to the case of repulsive interaction (U > 
0). Here we assume that h > and the energy dispersion 
is given by £a=i(2)(fc) = —2t cosk— {+)h — fx, where a = 1 
(a = 2) corresponds to the majority (minority) spin. 
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FIG. 8: (Color online) The spectral functions yl_R.i(g,aj) and 
AR,2{q,i^) for U/t = 3 and (J = 0.01/a, with fixed h/t = 0.2 
(top) and 0.6 (bottom). 



We denote the filling factor by p (0 < p < 2). For 
h < t(l — cosTrp), two bands are overlapped at Fermi 
energy. The Fermi momenta and the bare velocities are 
given by 




(4.44) 



(4.45) 



The effective velocities u± are given by Eq. (|2.2ip with 
g = U and Ki = K2 = 1. By considering incommen- 
surate band filling, the umklapp scattering is irrelevant. 
For /i = 0, the Hubbard Hamiltonian (I4.43P possesses 
the SU(2) symmetry and the analysis given in Sec. IIVAI 
applies. The neglect of the backward scattering can be 
justified for strong magnetic field, and the effective model 
can be described by Eq. (|2.ip with g = U. The magnetic 
field dependence of the bare velocities ui^2 and the renor- 
malized velocities u± are shown in Fig. [T] For h = 0, 
the velocities u+ and u_ correspond to the conventional 
charge Vp and spin v„ excitation velocities, respectively. 
The spectral functions Aii^a{q,oj) for h/t = 0.2 and 0.6 
are shown in Fig. [5] In Fig. [51 the contour plot oiAa{k, uj) 
is shown in the range —0.05 < (fc — kp^a) < 0.05 and 
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FIG. 9: (Color online) Contour plot of the spectral function Aa(fc,aj) for a = 2 (left) and a = 1 (right), with fixed U jt — 3 
and h/t = 0.2 (top) and h/t = 0.6 (bottom). The dotted lines denote the bare dispersion ea{k) ~ Ua{k — kp.a) where the Fermi 
momenta are given by {kF,2,kF,i) — (0.64,0.93) for h/t = 0.2 and (0.35,1.22) for h/t — 0.6. In the insets, the bare energy 
dispersions are shown and the rectangular represent the region where the spectral functions are plotted in the main figures. 



-0.05 < uj/t < 0.05, with fixed h/t = 0.2 (top figures) 
and 0.6 (bottom figures). The momentum k is related to 
qhy k = kF,a+q [see Eq. (|4.3p ]. For weak magnetic field, 
the spectral functions exhibit similar behavior obtained 
in the SU(2) symmetric case,-'*-'-- except for the non-zero 
weight seen at ~u^q < uj < —u^q. In the SU(2) sym- 
metric case, there is no weight at —u+q < uj < —U-q (see 
Fig. [3]). For strong magnetic field, the difference between 
Ml and U2 becomes large and the spin-up and spin-down 
electrons are effectively decoupled, where Ar^i {Ar^2) has 
stronger weight at oj w u^q {uj w u^q). We note that the 
weight of Aii^2{qi'^) at a; ~ u^q becomes relatively large 
compared with that of Aii^i{q,uj) at w ~ "+9, reflecting 
the large density of states for spin-down particle. Espe- 
cially, Afl^2(9,w) exhibits a similar behavior obtained in 
the spinless case [Fig. 1(a) in Ref . |49| . We also note that, 
as seen from Fig. [SI the two branch feature is prominent 
near Fermi energy (w ~ 0), however, one-particle feature 
is recovered in the high-energy region. As seen in Fig. 
[71 the velocity u+ takes a close value to ui, while u_ 
is strongly renormalized from the bare velocity U2. For 
the spectral function yl2(fc,w), a strong singularity can 
be seen &i lj — u^q, while relatively large weight can be 
obtained at — ju+gj < uj < [u+qj for small uj. 




By using Eqs. (|4.3ip and (j4.32p . the h dependence of 
the exponents /3s are obtained (Fig. ITUl) . As mentioned 



FIG. 10: (Color online) Exponents /?a,± and /3a,± as a func- 
tion of the magnetic field for U/t = 3 and p = 1/2. 



before, the spectral function has two peaks at uj — u-^^q 
and UJ = u^q for /?2,+ < and /?i,- < 0, respectively, 
while each divergence is replaced by a cusp for /32,+ > 
and /3i,_ > 0, respectively. Figure [TT] shows such a 
critical value of h/t as a function of U/t, for /?2,+ = 
(sohd line) and /3i_ = (dotted hue). The dash line 
corresponds to the upper bound of h/t for the stable u_. 
With increasing U/t, the critical value of h/t for /32,+ = 
and /3i._ = decreases. 
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ing function will be such that J"i/(0) is a constant, while 
Ti,{x) ~ a;'^"'" when x — >■ oo. A similar, but weaker ef- 
fect should be seen on the cusps at w = —u±q. Again, 
the power laws will be replaced by scaling functions that 
reproduce the infinite slope for T = 0. Finally, we ex- 
pect that the gap —u^q < u) < u^q will start to fill. For 
higher temperatures, T ^ |u+ — U-\q, we expect that 
the difference between the peaks at a; = u±q becomes 
blurred, and only a single broad peak will be observed. 
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FIG. 11: The f/-dependence of the critical value for h/t sat- 
isfying 132, + — (solid line) and /3i,_ — (dotted line), for 
p = 1/2. The dash Une denotes the upper bound of h/t for 
the real U2, i.e., corresponding to U2c- 



V. CONCLUSION 

In the present paper, we have derived expressions of the 
fermion spectral functions of a two-component Luttinger 
liquid at zero temperature in terms of Appell hypergeo- 
metric functions. We have shown that in the SU(2) sym- 
metric case, these expressions reduce to the Gauss hyper- 
geometric functions. Our expressions allow the recovery 
of the singularities derived in Refs. i49:.50 but also describe 
the amplitude of the singularities as well as the behavior 
of the spectral function away from the singularities. The 
results of the present paper could be used to calculate 
zero temperature fermion spectral functions of various in- 
tegrable models^"— using exact results on the TL-liquid 
exponents. In the case of three-component models,®'^ a 
more general version of the Feynman identity— is appli- 
cable, and we expect that the zero temperature spectral 
function will also be expressible in terms of generalized 
hypergeometric functions of three variables. In the case 
of non-zero temperature, the real-space Green's function 
is expressed as a product of powers of hyperbolic sines^ 
and the Feynman identity^'^ we used to derive our ex- 
pressions does not lead to a tractable expression. How- 
ever, it is still possible to discuss the qualitative changes 
to the spectral functions produced by positive temper- 
ature. The first effect of being at finite temperature is 
that the divergences of the spectral functions at lj — u^q 
{f = ±) are cut-off by the thermal length which is of 
the order of u^/T. Therefore, the peaks of the spec- 
tral function ^jj_a(g,a;) will be replaced by maxima of 
height ~ T^"^'" as w — )■ u^q- For \u) — u^q\ '^ T, we 
expect that the integral giving the spectral function is 
not strongly affected by the finite temperature as the 
largest contribution comes from integration over length- 
scales smaller than the thermal length. Therefore, for 
the lowest (non-zero) temperature scale, we expect that 
the spectral function will be modified by the replace- 
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Appendix A: Integrals 

1. Fourier transforms 

To obtain the spectral functions of the fermions, we 
need the Fourier transform: 



Jiq,Lu) = 



dxdt e-K92;-c^t) 



[a + i{uit - x)p^ [a + i{u2t + x)p-' 
in the linrit of a — > 0+. By the change of variables: 



,(A1) 



t = 



X1 + X2 

Wl + M2 



U1X2 — U2X1 
U1+U2 



(A2) 



we find: 



J{q,uj) 



1 



Ui +U2 



.■('i"2+") 



Xi 



dX, 



{a + iXip^ 



dXo 






(a + 1X2)'^^ 



(A3) 



The two integrals in the product are obtained from the 
formula (3.382.7) of Ref. 84, giving the final result: 



J{q. 



{2ttY {uj + U2qV'~^{uJ-Uiq)^-'-^ 
r(7l)r(72) {ui + U2V^+'^-'-'^ 

xe{uj-uiq)Q{oj + U2q), (A4) 



T^"" J"^[(a; - Ui,q)/T] in Eqs. (|i^ - (|i^ . The seal- where the hmit a ^ 0+ has been taken. 



2. Integrals expressible as hypergeometric 
functions 



The integral representation of the Gauss hypergeomet- 
ric function is 



X2Fi{a + (3 + j~l, (3;a + (3; 



b{- 


-a — 


c)\ 


a{' 


'b + 


(A8) 



2^1 (a, /3; 7; z) = Y(B)T( - B) / ^^ (1 - tz)°' ' where we have used Eqs. (K5\ and (|M)) . 

(A5) 
which satisfies [Eq. (15.3.3) in Ref. ^ For the evaluation of Eq. (I4.35p . we meet the integral: 

2F,{a, /3; 7; ^) = (1 - z)^-"-^Fi(7 - «, 7 - /3; 7; ^). 

(A6) 
Typical integral formula for obtaining the spectral „i 

functions Km is given by ^1 = / '^^ (1 - tf'H^'-\Fi{a, ^ + ^'; c; A + a**)- (A9) 

Jo 



^-r"--"-)'- ,^:;::^:::-. - <-> 



By changing variable t = 1 - c(l - w)/(c - &w), we find By changing the variable i by t = Xi^/{X + fJ. - /iz^), and 

(a + 5)/3+7-i r(a)r(/3) 
a'3(-6 + c)°+'3+'r-i r(a + /?) 



(,+,)...-! r(a)r(/3) i"^, ^^p;;f ::^ ,f;/«' '■' ^^ !) r ^^- "^^ "'"^ 

-^ = „fl^ u I „N»+/9+-v-i T^^„ , p\ (a)„ = r(a + n)/r(a), wefind 



,f^o^ ^ (c)„ «! io (A + Ai - H''+'' +" 

_ r(/3)r(/?') / A \'''^ 1 (a)„(/3 + /3')„ ^ /« , «, , «,.« , «,. M 






r(/? + /3') VA + /i; ,^„n! (c) 

r(/3)r(/3') ^ A N^^' ^ (/3 + /3')„+™(a)„(/3')mA"(Af^r 



r(/? + /3') VA + /i; „^o (c)„(/3' + /3U n\ m\ 

where i^2 is the Appell hypergeometric function. By using the relation: F2{A; B, B';C, A;x,y) = (1 — 
y)'^'Fi {B; A - B', B'; C; x, x/{l - y)), we obtain: 

h = B(/3,/3')i^i(a;/3,/3';c;A,A + A*), (All) 

where Fi is the Appell hypergeometric function. Similarly, the integral in Eq. (|4.33p can be performed as: 

h = f dt{l-t)'^-H'^'-^(l + ^ty \fi (/? + /?' + c-l,&;c;A + Mi) 



vf^ (^±i^:±4^iM^^ r..(i-.)^-(.)^'-fi-^ 

;^VA + M/ (c)n n! V A + ^7 



A \ n^o o'\ Y^ A" (a+^) ,„ . of , i\ (^)n(/3') 

B[fi,(i) > — j- ^(/^ + /3 +c-l)„^ 



l + ?7l ~ 



A + /.; -'-'^Z^^n! m! ^ ' '^ ' ^"-'" (c)„(/3 + /3')™ 

-^) S(/3,/3')i^2f/3 + /3' + c-l;6,/3';c,/3 + /3';A,-^y (A12) 

A + /i/ V A + /i/ 
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Appendix B: Finite size bosonization and exponents momentum operator, 



It is well known in conformal field theory that there is 
a relation between the energy-momentum tensor and the 
Virasoro generators that give the conformal weights (or 
dimensions) of the operator a^^i^^ . As a result, it is possi- 
ble to relate the dimension of a given operator to the en- 
ergy of the state generated by acting on the ground state 
of the Hamiltonian of a conformally invariant model with 
that operator—. Multicomponent models are not in gen- 
eral conformally invariant, but their critical properties 
can be obtained by considering a semi-direct product of 
Virasoro algebras^. A more elementary approach, that 
we will follow here, uses finite size bosonization^*'^^. In 
that approach, the fields <j)a and 9a in Eq. (|2.4p admit 
the decomposition: 

M^) = 0^"^-^a;^--^^(/.,(g)e^'^^ (Bla) 

q 

daix) = e^^'>-^x + ^Y.^a{q)e'''^. (Bib) 

In Eqs. (JB1[) . we have the commutation relations: 
[i?!)q , Jb] — —iSab and [Oq\ rib] — —i^ab- In the Hamilto- 
nian, (pQ and 9q do not appear as a result of the deriva- 
tions, but there is an extra term: 



2i 



SH = —y^^MabJaJb + NabngUb) 

*nNn), (B2) 



2i<'-"" 



[J = *(Ji, J2) and n — *(rii,ri2)] which vanishes when 
L — > 00 and is called the zero mode contribution. The 
ground state has Ja = and n^ = 0. If we consider the 

I 



-P == ^ / dxllad^(l)a, 



(B3) 



it also contains a zero mode contribution equal to: 

SP^jY. naJa - j'nJ. (B4) 



If we use the relations M = Q{A2)^^'^*Q and N 
P(A2)^'^^*P, we can rewrite: 



SH = — 



(*JQ + *nP)(A2)i/2(*QJ + *Pn) 
(*JQ-*nP)(A2)i/2(*QJ-*Pn)l , (B5) 



6P = [(tjQ + tnP)CQJ + 'Pn) 

-CJQ-'nP)CQJ-*Pn)] . 



(B6) 



Now, if we act on the ground state with the opera- 
tor e*^"^''''^""'"^''''^'')^'^^ due to the presence of the term 
J2a(''la^o + Ca'/'o )' ^^^ resulting state will belong to 
the subspace with tIq = —rjajJa = — Ca- Its zero-mode 
contribution to the ground state energy will be: 



SP = ^[C^Q + *vP)l-C^Q-'vP)l] 



(B7) 
(B8) 



For an infinite system, a simple generalization of Eq. 
shows that the correlation function take the form: 



/e*(l^a'?-«o+€a0a)(a:,t)g- 



»(E„';ae„+?„</.j(o,o)^ ^ tt 



a 



a + i{ui3t + x) 



CvP+*iQ)l/4 




CvP-iQ)l/4 



(B9) 



So we see that the dimensions {^rjP ± *^(5)o/4 in the 
correlation functions also appear in the zero mode con- 
tributions to the excited state energy and momentum of 
the finite size system. This remark is the basis for the 
method of Frahm and Korepin>2^ Indeed, Eq. (8.2) of 
Ref. [8^ is recovered with P = *C/*(Z~^) where: 



U = 



1 1 
1 



(BIO) 



since Frahm and Korepin have defined Nc = N-^ 



Ni 



and Ns = N^, so that: 



N, 






(Bll) 



Similarly, Eq. (3.6) of Ref. 188| is recovered by taking 
Q ~ U^^Z. Using these definitioiis, we can recover Eqs. 
(3.12), (3.11) and (3.13) of Ref. m These relations 
have also been derived in Ref. |72 by comparing the crit- 
ical exponents derived from the Bethe Ansatz with the 
ones derived from bosonization. 
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